Creation of vortices in a Bose-Einstein condensate by a Raman technique 
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We propose a method for taking a Bose-Einstein conden- 
sate in the ground trap state simultaneously to a different 
atomic hyperfine state and to a vortex trap state. This can 
be accomplished through a Raman scheme in which one of the 
two copropagating laser beams has a higher-order Laguerre- 
Gaussian mode profile. Coefficients relating the beam waist, 
pulse area, and trap potentials for a complete transfer to the 
m = 1 vortex are calculated for a condensate in the non- 
interacting and strongly interacting regimes. 
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Based on the rapid experimental progress in the field 
of atomic Bose-Einstein condensation jl]-^ , the prospects 
look good that experimenters will soon be generating and 
observing vortices in these sytems. Vortex solutions to 
the Gross-Pitaevski equation for the condensate wave- 
function have already been found numerically iQjHj. It 
has been pointed out that unlike homogeneous Bose flu- 
ids such as ^He, there exist two distinct concepts of sta- 
bility in a trapped Bose-Einstein condensate (BEG) [^. 
Although they are thermodynamically stable while the 
system is being rotated, this does not necessarily imply 
dynamical stability once the rotation ceases. Thus, one 
might be interested in the dynamics of the vortex. A 
necessary step is to start with a well-defined state; here 
we propose an optical method for coherently creating a 
vortex state from the ground state. 

It is well known that the higher order Laguerre- 
Gaussian modes of a laser resonator possess orbital angu- 
lar momentum . (Note that this is possibly in addition 
to the angular momentum contributed by circular polar- 
ization of the light) . If this angular momentum can be co- 
herently transferred to a BEG, the circular motion of the 
atoms should result in a vortex state. Since this involves 
removing photons from the Laguerre-Gaussian mode, we 
need to avoid subsequent spontaneous emission and heat- 
ing of the condensate 1^. One way to accomplish this is 
to use a Raman scheme in which the second beam is an 
ordinary Gaussian mode, and the atoms make a transi- 
tion between hyperfine states along with the transfer of 
angular momentum. To transfer all of the atoms from 
the trap ground state to the vortex state, we need to find 
the equivalent of a two-photon tt pulse. 

We denote the condensate wavefunctions for hyperfine 
states 1 and 2 by ipi and "02. The trapping in current ex- 
periments is done magnetically, so we allow for different 
trap frequencies loi, loi for the two atomic states. Since 



the vortex can already be described in two dimensions, 
we will treat the condensate as though it simply has a 
thickness L in the z-direction. The light beams both en- 
ter from the top (or bottom) and are copropagating. An 
important assumption is that the condensate be optically 
thin in this dimension so that the light beams will affect 
all of the atoms. For the time being we neglect the colli- 
sions between atoms; their effect will be explained below. 
It is convenient to define 
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Then, after adiabatically eliminating the third level from 
the Raman transition [Bj, the modified Gross-Pitaevski 
equations are given in polar coordinates by 
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Note that we use unit normalization. 
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The two-photon Rabi frequency is 
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and in the unsealed variables 

f^(x,0=(^^-^^("'g2''''''^"''^ (12) 

where di,d2 are the dipole matrix elements between 
states 1 and 3, 2 and 3 respectively, and Ei,E2 are the 
near-resonant fields detuned by A from those transitions. 
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The fields have oppositely circular polarizations as ap- 
propriate to a miT- to mp^2 transition. We will choose 
the field Ei to be an ordinary Gaussian mode (TEMoo), 
while E2 will be the TEMn Laguerre-Gaussian mode. 
(In fact it makes no difference which of the two fields has 
which mode structure.) We approximate these modes by 
their values at 2 = 0, 



Ei(x,t) = Ai(t)exp [~r^/wl] 
E2 (x, t) = A2 {t) — exp [~r^ /w: 

W2 



Substituting into Eq. (|l^) yields 
fi(x, i) = r2o(i) — exp 

where 
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always be unoccupied. Second, when collisions between 
the atoms are included (see below), the differences in 
chemical potentials for higher states will no longer be 
on resonance so the transitions are weakly driven in any 
case. Substituting into the Gross-Pitaevski equations (|^) 

and using the orthogonality of the V'j™'' ; we find the equa- 
tions of motion 
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where 
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and 



(di-Ai(t))*d2-A2(0 
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These are exactly equivalent to the Bloch equations for 
a two-level system. We note that an additional detuning 
/ij — Ml"'' is required from the untrapped atom reso- 
nanant Raman transition; the energy this represents is 
imparted to the gas to produce the rotation. The effec- 
tive Rabi frequency is given by 



We begin with all the atoms in hyperfine state 1 and 
the associated trap ground state. The ground trap states 
of each hyperfine state are (we drop the bars on the scaled 
variables) 
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and the two vortex states rotating clockwise and coun- 
terclockwise are 



g{t) = noit)La^ 
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As usual ||T^, one may obtain a complete transfer to the 
upper level at time i/ by choosing an on-resonance tt- 
pulse. 
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We can expect that the technique will remain ap- 
plicable when interactions between the atoms are in- 
cluded. We replace the ground and vortex wavefunctions 
in Eq. (|2^) by their Thomas-Fermi approximations in 
the strongly interacting regime. 



We expand the time-dependent condensate wavefunc- 
tions in this basis. 



V',(x,t)=cf)(t)V.;"^(x)e 
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for J = 1,2 
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where 
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Here each ^j™'' is the chemical potential when all atoms 
are in the hyperfine state j and the condensate has az- 
imuthal angular momentum m. We have neglected the 
higher terms of the series for two reasons. The Raman 
transition couples these higher trap states which are ini- 
tially empty only to each other and therefore they will 



(32) 



is the scaled interaction potential between two atoms in 
hyperfine state j. These functions then appear in the 
integral for the effective Rabi frequency Eq. (pa). (Note 
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that for the vortex, the Thomas- Fermi approximation in- 
cludes the azimuthal kinetic energy but neglects the ra- 
dial part.) The chemical potential of the ground state is 
found from normalization to be 
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To find the chemical potential of the vortex state, we 
must solve numerically the normalization equation 
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and then set 



^2"^^ — coth(^. 
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It is desirable to make the effective Rabi frequency 
large, so that short pulses can be used and the dynam- 
ics of the condensate can be ignored while the vortex is 
being made. To this end, we numerically evaluate the 
integral Eq. (|2^) and determine its dependence on the 
waist wi and W2 of each beam. We choose parameters 
similar to those of the JILA experiment with overlap- 
ping condensates in two hyperfine states For a total 
of = 2 X 10® condensed ^^Rb atoms, mean of trap 
frequencies ^^1^2 = 2tt x 400Hz, thickness L = lO^^m, 
a — -\/2, and scattering lengths ajj — lO^^m, we ob- 
tain the ground and vortex states shown in Fig. |l|. (The 
mean trap width is a± = 540nm.) The required detuning 

{iJ,^^ —^[~^^)/h = 39.4kHz is mostly due to the difference 
in trap potentials in this example. In Fig. ^jwe show the 
dependence of g/rig on the beam waists. We note that 
there is a maximum as a function of W2 when wi is held 
fixed, for a value of W2 on the order several a±. It is 
advantageous to choose wi as large as possible; for the 
case of a plane wave Ai, the optimal waist size for the 
Laguerre-Gaussian beam is 1^2 = 19.3aj^. 

The resulting vortex state has a zero condensate den- 
sity at its center. However, current techniques are prob- 
ably not capable of imaging this small core in the atomic 
density directly. In liquid ^He, a negative ion can be 
trapped at the vortex core, and after accelerating the ion 
out of the fluid it can be imaged on a phosphor screen 
PH . It might be possible to adapt this technique to the 
atomic BEG. An alternative way of observing the vor- 
tex is to allow it to interfere with another (non-vortex) 
condensate, and image the fringes as in A similar 
technique has been used in nonlinear optics [T^ . The 
27r phase shift around the vortex will show up as two 
interference fringes merging into one. 

We have shown how a vortex state can be created in 
a trapped Bose-Einstein condensate. The transfer to the 
new state is completely coherent, as it involves only stim- 
ulated optical transitions. It could also be applied to 
create higher vortex states with |m| > 1. To investigate 
persistent currents, one may transfer the atoms to the 
vortex while a repulsive potential is placed at the center; 



a simple way to do this is to detune the Gaussian laser 
mode above the single-photon resonance and leave that 
field on after the Laguerre-Gaussian pulse has ended. 
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FIG. 1. Thomas-Fermi wavefunctions of the initial and fi- 
nal condensate, appropriate for strongly interacting atoms, 
for the parameters given in text. Solid curve: Ground trap 
state, hyperfine state 1 Dashed curve: Vortex trap state, hy- 
perfine state 2. Note this function has a small core at which 
the density goes to zero. 
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FIG. 2. Overlap integral g/0.o as a function of La- 
guerre- Gaussian beam waist W2 for wi = oo (solid curve), 
wi = 50ox (long dash) and w\ = 20ox (short dash). 
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